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The Letter considers the ground state and the Tkachenko modes for a rapidly rotating Bose- 
Einstein condensate (BEC), when its macroscopic wave function is a coherent superposition of 
states analogous to the lowest Landau levels of a charge in a magnetic field. As well as in type II 
superconductors close to the critical magnetic field H C 2, this corresponds to a periodic vortex lattice. 
The exact value of the shear elastic modulus of the vortex lattice, which was known from the old 
works on type II superconductors, essentially exceeds the values calculated recently for BEC. This 
is important for comparison with observation of the Tkachenko mode in the rapidly rotating BEC. 
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A rapidly rotating Bosc-Einstein condensate (BEC) of 
cold atoms is now a subject of intensive experimental and 
theoretical investigations HSU SUB 01111 It is well 
known that rotation gives rise to a regular triangular vor- 
tex lattice. At moderate rotation speed this is a lattice of 
vortex lines (points in the 2D case) with the core size of 
the order of coherence length £, which essentially smaller 



than the intervortex distance b — y n/y/Sft. Here is 
the angular velocity of rotation and k = h/m is the circu- 
lation quantum. One can call it the Vortex Line Lattice 
(VLL) regime. With increasing fi, the vortex lattice be- 
comes more and more dense and eventually enters the 
regime, in which the vortex cores start to overlap, i.e., 
£ becomes larger than the b. This regime is analogous 
to the mixed state of a type-II superconductor close to 
the second critical magnetic field H C 2 ~ &o/£, 2 ($o is the 
magnetic flux quantum), at which the transition to the 
normal state takes place. However, in a rotating BEC 
there is no phase transition at corresponding "critical" 
angular velocity il C 2 ~ k/£ 2 . Instead the crossover to 
the new regime takes place: At f2 3> f2 C 2 ah atoms con- 
densate in a state, which is a coherent superposition of 
single-particle states similar to the Lowest Landau Levels 
(LLL) of a charged particle in a magnetic field (the LLL 
regime). An important method of investigation of the 
vortex structure is studying its collective modes. Cod- 
dington et al. [jj] were able to detect the Tkachenko 
modes (transverse sound in the vortex lattice) in the VLL 
regime experimentally. Recently they increased the ro- 
tation speed in the attempt to reach the LLL regime . 
They revealed softer Tkachenko modes as was predicted 
by the theory for the LLL regime 



Theoretical study of the Tkachenko mode requires 
good knowledge of the equilibrium state. A number of 
papers addressed this issue using the analogy with the 
quantum Hall effect H B H H @ . They started from 
the LLL wave functions for noninteracting particles in 
a trapping potential and switched interaction on after 
it. For a regular vortex lattice this yielded the Gaus- 



sian density profile |5| but it was unstable with respect 
to small distortions of the lattice, which transformed the 
Gaussian profile to the inverted parabola (Thomas-Fermi 
distribution) BBBBEl This Letter suggests another 
strategy. One can start from an infinite periodical vor- 
tex lattice in an infinite uniform liquid, neglecting first 
the trapping potential but taking into account interac- 
tion. The exact wave function for this state was found in 
the classical work by Abrikosov for type II supercon- 
ductors close to H C 2, and later it was gen eralized for an 
arbitrary unit cell of the vortex lattice [111 |l2j . If a trap- 
ping potential is added, the vortex lattice (as well as the 
liquid density) ceases to be uniform. But as far as mod- 
ulation by the trapping potential is smooth (the cloud 
size is much larger than the intervortex distance) it does 
not distort the lattice essentially (except for the extreme 
periphery of the cloud) and one can use the thermody- 
namic potential derived for a uniform vortex lattice. The 
suggested approach is especially useful for investigation 
of the Tkachenko mode since for the infinite uniform lat- 
tice the shear elastic modulus can be calculated exactly 
and for the vortex lattice in type II superconductors it 
was done many years ago |l3| . The previous calculations 
of the shear modulus in BEC B m yielded the values 
smaller than the exact one. This difference is important 
for comparison with recent experiments on Tkachenko 
modes |2(. . 

We consider a 2D rotating BEC in a parabolic trap- 
ping potential characterized by the frequency uj± . In the 
Gross-Pitaevskii theory the Gibbs thermodynamic po- 
tential is 
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Here tf> is the BEC wave function, fj, is the chemical po- 
tential, g is the interaction constant, and vo = [H X r\ 
is the velocity of solid body rotation. The Gibbs poten- 
tial is invariant with respect to the gauge transformation 



if) — > -0 e j — * + (k/27t)V0. At rapid rotation the 
potential mil 2 r 2 /2 of centrifugal forces nearly compen- 
sates the trapping potential muj 2 _r 2 /2, Though stability 
requires that oj± > Q, at the first stage of the analy- 
sis one can assume that lj±_ = Q and BEC is infinite in 
size. Then the Gibbs potential Eq. is invariant with 
respect to translation accompanied by the gauge trans- 
formation, which corresponds to the shift of the rotation 
axis. 

As well as for type-II superconductors close to H c2 , in 
zero-order approximation one can neglect interaction in 
the LLL regime £ 3> b. Then the linearized Schrodinger 
equation is similar to that for a charged particle in a 
uniform magnetic field: 
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At fi — hfl it has a solution, which corresponds to the 
lowest Landau level: 
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where I 2 — k/4ttH and yk = —l 2 k. The solution is given 
for the gauge with VQ(—2Qy, 0). The frequency 2f2 is the 
analog of the cyclotron frequency tu c — eH/mc for an 
electron in a magnetic field. If we consider a square Lx L 
with periodic boundary conditions, then k = —2im/L 
with the integer n. Using the condition < yk < L, one 
can see that the integer n should vary from zero to the in- 
teger closest to L 2 /2nl 2 . This is the total number of LLL 
states, and the density of the LLL states is l/2nl 2 . All 
these states are orthogonal each other and have the same 
energy. But degeneracy is lifted by taking into account 
the interaction energy. The solution, which corresponds 
to the periodic vortex lattice with one quantum per lat- 
tice unit cell, is Il2j 
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where C„+i = C n exp(2n ib cos a / a) , a, 6, and the angle 
a are the parameters of the unit lattice cell (see Fig. 

This solution yields the thermodynamic potential of 
the infinite BEC in the LLL regime: 
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where n = M\ 2 ) is the average particle density and the 
parameter |l2j 
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characterizes dependence on lattice parameters a, b, 
a via the complex parameter 
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FIG. 1: Vortex lattice before (solid lines) and after (dashed 
lines) shear deformation 



Here 
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are theta functions [lij . The minimum of the interaction 
energy correspond to the triangular vortex lattice with 

[3 = 1.1596, a = b = 21^JttJ\^, a = n/3. According 
to Eq. JSJ the Gibbs potential has a minimum at the 
particle density n = {ft — hfl)/ ' [3g. The vortex density 
n v — \/2tt1 2 is equal to the density of LLL states. 

Let us now take into account the parabolic trapping 
potential. In equilibrium the usual Thomas-Fermi con- 
dition takes place: 
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where fi(0) is the chemical potential in the center of the 
trap, r = 0. This gives the inverted parabola for the 
density distribution in the trap: 



n(r) = n(0) 

The condition n(R) 
c s {0)VW(u± 
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yields the cloud radius R 



fl 2 ), where c s {0) = ^/3gn{0) is the 
sound velocity in the center (r = 0) of the rotating BEC. 

The authors of Refs. 0,0 also received the inverted- 
parabola density profile. In their approach it looks as 
the effect of small distortions of the vortex lattice, which 
essentially transforms the Gaussian density profile pre- 
dicted by Ho [5j for an ideally regular vortex lattice in 
a trap. But sensitivity of the Gaussian profile to small 
distortions is nothing more than evidence that the Gaus- 
sian profile is unstable. Small distortions are inevitable 
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in a finite vortex cluster (see Sec. V.B in Ref. 0), but 
it does not mean that they are crucial for bulk proper- 
ties of a macroscopic vortex crystal. Our approach yields 
the inverted-parabola density profile without paying at- 
tention to small distortions, because the Gaussian profile 
does not appear at any stage of the analysis. 

Let us consider restrictions on existence of the LLL 
regime. First, the energy of the lowest Landau level, 
hfl, should exceed the interaction energy per particle, 
(3gn (but is close to unity and is not essential for an 
order-of-magnitude estimation). This yields the inequal- 
ity n <C fifl/g, which is equivalent to the inequalities 
£ > b or ct > c s 



where c s = y/gn and ct = y k£lI%ti 
are the sound and the Tkachenko mode velocity in the 
VLL regime. Second, the BEC with a regular vortex lat- 
tice exists as far as the filling factor n/n v (the number 
of particles per vortex) exceeds unity (see below), i.e., 
the inequality n 3> 0,/k is required. The two inequali- 
ties are compatible for a weakly interacting bose gas with 
g -C h 2 jra. Since g ~ h 2 a s /ml z , it is necessary that the 
oscillator length l z for the trapping potential localizing 
the BEC cloud along the rotation axis exceeds the scat- 
tering length a s . One can rewrite these inequalities in 
terms of the total number of particles TV = unR 2 . Since 
R 2 = 2(3gn/m(uj 2 _ — fi 2 ), the order-of-magnitude estima- 
tion for density is n ~ 
regime takes place if 



Eq. (|13|l . one obtains the value of the shear modulus: 
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One more condition is the presence of many vortices in 
the cloud: im v R 2 3> 1. This yields the inequality 
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gm 
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This is compatible with the previous inequalities only for 
rapid rotation when lo\ — £1 2 -C SI 2 . 

Deformation of the triangular lattice should increase 
its energy. On the basis of the elasticity theory for a 2D 
crystal with hexagonal symmetry , the elastic energy 
should be 
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■l = C\(U XX + Uyy) 2 + C 2 [{U XX ~ Uyyf + 4tt^ y ] (13) 



where C\ is the compressional modulus, C2 = mncip/2 
is the shear modulus of the vortex lattice, Ct is the 
Tkachenko wave velocity in the LLL regime, and tiy = 
\{y%Uj + V jUi) are components of the deformation ten- 
sor. Since the parameter (3 does not depend on the vortex 
density, the compressional modulus vanishes: C\ = 
In order to find the shear modulus let us deform the tri- 
angular lattice as shown in Fig. ^ The variation of the 
complex parameter £ is proportional to the shear defor- 
mation: <5£ = Sp = 2u xy sin a. Expanding the expression 
Eq. © for j3 and comparing the term cx Sp 2 to the ther- 
modynamic potential, Eq. (jSJl, with the elastic energy 
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This exactly agrees with the numerous calculations of the 
shear modulus cqq = 2C2 for the flux lattice in type II 
superconductors close to the critical field H C 2 0] . But in 
this work it was calculated anew since Eq. (|14|l yields the 
result 2 times larger than the value by Sinova et al |3J and 
10 times larger than the value (81/807r 4 )rrmc 2 received 
by Baym 0. This numerical difference is important for 
interpretation of the experiment (see below) . 

One can proceed with the analysis of the Tkachenko- 
mode spectrum on the LLL regime using the same hy- 
drodynamic equations as in the VLL regime [ll[ll[l3, 
but with c s and Ct replaced by c s and ct- However, 
since the LLL regime is possible only for very rapid rota- 
tion when oj± — fl -c oj± the effect of high BEC com- 
pressibility is always strong (see the comment [16] in 
Ref. 0) transforming the Tkachenko mode spectrum to 
quadratic: ui — c s CTk 2 /2Q. In this limit one can use 
the analytic expression for the Tkachenko eigenmodes 
of a finite BEC cloud in a parabolic trap derived for 
the VLL regime in Ref. [l8f &i = 7</s. Here tOi is 



the value of the reduced frequency Cj — y \/3 / 8ttlu R / ct 

and s — 2\/2Q/ y 1 uJ 2 ^ — Q 2 . The numbers ji depend on 
the number i of the eigenmode. The two lowest eigen- 
modes correspond to 71 = 7.17 and 72 = 16.9. Using 
the expression R = y/2c s (0)/y/u)\_ — ^ 2 for the cloud ra- 
dius one obtains uJi — \J it /y/^^i (ct/c s )(o; 2 _ — fi 2 )/fi 



2W tt I V37i(cT /c s )(^j_ — fi). Applying this expression 
to the LLL regime one should replace c s and ct with 
Cs = y/j3c s and ct given by Eq. ifTfjl . This yields the 
Tkachenko eigenfrequencies w 0. 86337; (wj^ — O). The 
experiments on rapid rotation roughly agree with ex- 
pected linear dependence of the first Tkachenko mode on 
lu± — Q, but the slope of the experimental line is about 
4 times less than our prediction for the LLL regime. It 
could be an evidence that the experiment has not yet 
reached the LLL limit. Since experimental values of 
(wj_ — f2)/a;j_ look small enough, apparently in order to 
approach to the LLL limit further, the experiment should 
be done with a smaller number of atoms. 

What should happens with the LLL regime when the 
filling factor n/n v approaches unity? Possible answers to 
this question are now intensively studied by theoreticians 
numerically and analytically 0,0, El- The following sim- 
ple discussion certainly cannot be a substitute of these 
investigations, but hopefully would be useful for better 
qualitative understanding of possible scenarios. At low 
n/n v bose condensation is destroyed and a single macro- 
scopic wave function cannot describe all atoms anymore 
Q . For an illustration let us construct a state of the bose 
system, which successfully competes with the BEC state 
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at low filling factors. The single-particle state 
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corresponds to the classical Larmor orbit of a charged 
particle centered in the point with coordinates (Xi,Yi). 
It is a solution of Eq. J5J) for the gauge vo(— fly, fix) and 
belongs to the space of the LLL states. Let us assume 
that any LLL state given by Eq. (|15|l cannot have more 
than one atom (the number N of them is much less than 
the number of LLL states). Then the wave function of 
N bosons can be written as 



1 N 



(16) 



p i=i 



Here i is the number of the particle (subscript of x and y) 
and the number of the occupied site (JQ, Yi) (subscript of 
ip). In contrast to BEC, all atoms are in different states, 
and proper symmetrization (summation over all permu- 
tations P) should be done. Strictly speaking the states 
ipi given by Eq. (|15fl are not orthogonal. Therefore the 
wave function ^ requires some normalization factor. But 
for small filling factors overlapping of the states is weak 
and an additional normalization factor is exponentially 
close to unity The next step is to calculate the inter- 
action energy for this state, Only interaction between 
closest neighbors is essential. The cross interaction term 
between them is 
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Here ro is the distance between two sites. Assuming 
the triangular-lattice ordering of occupied sites, ro = 

^2/^/3n. Collecting all terms of this value we receive 
for the total interaction energy integrated over the whole 
space: 



6NM g , r 2,,2 3 



^■v/^/3n 



(18) 



Though we used the properly symmetrized wave func- 
tion of N bosons, in fact statistics is not important and 
the same energy can be derived for an unsymmetrized 
function. Comparing the energy Eq. (|18f) with the total 
energy of the BEC state, Eb = ^g(3Nn, one sees that at 
small filling factors n/n v <C 1 the BEC state has a larger 
energy. The transition between two states is expected 
at n/n v ~ 1, but the present calculation becomes inac- 
curate there. This estimation certainly cannot pretend 
to be an evidence that some crystal structure appears 
at small n/n v , but it demonstrates that in this limit the 
BEC is not the equilibrium state. 

Another option at n/n v ~ 1 is vortex melting with- 
out destruction of BE condensation 0, 0, 0. The 
plane Tkachenko mode u(k)e lkx ~ luJt , where u{k) is the 



amplitude of vortex displacements, has the energy ~ 
L 2 u(k) 2 mnc\k 2 in the square L x L. Considering the 
quantum melting at zero temperature this energy is equal 
to hu)(k)/2. The average vortex displacement squared is 
obtained by integration of u(k) 2 over the Brillouin zone 
of the vortex lattice: 



L z I kdku(k) 2 ^ / kdk X > 
/n mncTpk 1 



According to the Lindemann criterion, quantum melting 
is expected, when */ (u 2 ) approaches to the intervortex 
distance b ~ 1/ ^/n^. This corresponds to the filling fac- 
tor n/n v of order unity. So vortex-lattice melting can 
precede destruction of the BEC. 

In summary, the ground state of a rapidly rotating 
Bose-Einstein condensate is analyzed in the LLL regime 
using the exact wave function known for type II super- 
conductors close to the upper critical magnetic field H c2 . 
The analysis yields the inverted-parabola density distri- 
bution in a parabolic trap and the exact value of the shear 
elastic modulus of the vortex lattice in the LLL regime, 
which exceeds the values received in the previous calcu- 
lations for BEC. This has an impact on interpretation of 
recent experiments on very rapid rotating BEC. 

I thank Alexander Fetter for interesting stimulating 
discussions of the LLL regime. Ernest-Helmut Brandt 
helped a lot to my search of literature on vortex elasticity 
in type II superconductors. 
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